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Abstract 

We investigate finite energy solutions of the Einstein- Yang-Mills-Chern-Simons system in odd space- 
time dimensions, D = 2n + 1, with n > 1. Our configurations are static and spherically symmetric, 
approaching at infinity a Minkowski spacetime background. In contrast with the Abelian case, the con- 
tribution of the Chern-Simons term is nontrivial already in the static, spherically symmetric limit. Both 
globally regular, particle-like solutions and black holes are constructed numerically for several values of 
D. These solutions carry a nonzero electric charge and have finite mass. For globally regular solutions, 
the value of the electric charge is fixed by the Chern-Simons coupling constant. The black holes can 
be thought as non-linear superpositions of Reissner-Nordstrom and non-Abelian configurations. A sys- 
tematic discussion of the solutions is given for D — 5, in which case the Reissner-Nordstrom black hole 
becomes unstable and develops non-Abelian hair. We show that some of these non-Abelian configurations 
are stable under linear, spherically symmetric perturbations. A detailed discussion of an exact D = 5 
solution describing extremal black holes and solitons is also provided. 



1 Introduction 

In recent years the interest in the properties of gravity in more than D = 4 spacetime dimensions has 
increased significantly. This interest was enhanced by the development of string theory, which requires a 
ten-dimensional spacetime, to be consistent from a quantum theoretical viewpoint. Even in the absence of 
matter, solutions to the Einstein equations in dimensions higher than 3 + 1 exhibit properties which are 
strikingly different. For example, in a 4 + 1 dimensional asymptotically fiat vacuum spacetime with a given 
ADM mass and angular momentum, the geometry need not necessarily be that of the Myers-Perry [T| black 
hole. Notably, in this case there is the black ring [2] solution whose horizon topology is S 2 x S 1 , in contrast 
to S 3 of the former [T] . 

The rapid progress in the last decade has provided a rather extensive picture of the landscape of solu- 
tions for the five dimensional case [3], including configurations with Abelian matter fields [I]- [5]. Although 
the situation for D > 5 is more patchy, analytical [10] and numerical (TT| results suggest that the non- 
standard solutions found in D = 5 have higher dimensional generalisations; moreover, even more complex 
configurations are likely to exist as the spacetime dimension increases. 

In the case of higher dimensional gravitating systems of nonlinear matter fields, in particular with non- 
Abelian (nA) gauge fields, black hole and regular solutions are still relatively scarcely explored. This is an 
important direction since the theory of gravitating nA gauge fields can be regarded as the most natural 
generalisation of Einstein-Maxwell theory. Moreover, for the better known case of a D = 3 + 1 dimen- 
sional spacetime, the results in the literature show that various well-known, and rather intuitive, features 
of self-gravitating solutions with Maxwell fields are not shared by their counterparts with nA gauge fields. 
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For example, the Einstein- Yang-Mills (EYM) equations admit black hole solutions that are not uniquely 
characterised by their mass, angular momentum and YM charges, thus violating the no-hair conjecture |12| . 
Therefore the uniqueness theorem for electrovacuum black hole spacetimes ceases to apply for EYM systems. 
Also, in contrast with the Abelian situation, self-gravitating Yang-Mills (YM) helds can form particle-like 
configurations 13 . Another surprising result is the existence of nA solutions which are static but not spher- 
ically symmetric |14j . However, since it turns out that all these asymptotically flat solutions are unstable, 
their physical relevance is obscur43- (A detailed review of D = 4 gravitating particle-like and black hole 
solutions with nA gauge fields can be found in Ref. [15].) 

The study of D > 4 black hole solutions with non- Abelian matter fields is only in its beginnings. Based 
on the experience with Einstein-Maxwell solutions, it is natural to expect that higher dimensions D > 4 
allow for a rich landscape of solutions that do not have four dimensional counterparts. At the same time, 
considering such configurations is a legitimate task, since the gauged supersymmetric models generically 
contain non- Abelian fields. 

Most of the solutions displayed so far in the literature are spherically symmetric (an exception, being the 
results in |19j). As a new feature and in contrast with the situation in the D = 4 case, a generic property 
of the asymptotically flat higher dimensional EYM solutions is that their mass and action, as defined in 
the usual way, diverge [20] , [21], [21]) [22] • This can be understood heuristically by noting that the Derrick 
scaling requirement [24] is not fulfilled in spacetimes for dimension five and higher. Finite energy solutions 
exist only when the usual YM system is augmented with higher derivative corrections in the nA action 
[25) . Such terms can occur in the low energy effective action of string theory and represent the gauge field 
counterparts of the Lovelock gravitational hierarchy [26] (for a review of these aspects, see [27]). 

In a recent work [28j , a different way of regularising the mass of a D = 4 + 1 dimensional asymptotically 
flat, gravitating nA solutions was proposed. This was done by introducing a Chern-Simons (CS) term in the 
action. The CS density is a higher order term in the YM curvature and connection, and as such can be viewed 
as an alternative to the higher order curvature terms of the YM hierarchy employed previously in |22[ 125) . 
It turns out that this prescription does result in finite mass globally regular and black hole solutions and 
leads to a variety of new features as compared to the well known case of D = 4 EYM solutions. For example, 
these configurations cary an electric charge, emerging as perturbations of the Reissner-Nordstrom (RN) black 
holes. Moreover, in contrast to all other known asymptotically flat nA black holes without scalars, some of 
these solutions in |28) were found to be stable under linear, spherically symmetric perturbations. Also, for a 
particular value of the CS coupling constant, it was possible to construct both solitons and extremal black 
hole solutions, by exploiting the model in [29]. 

In this work we propose a general framework for the study of Einstein- Yang-Mills-Chcrn- Simons 
(EYMCS) solutions for an arbitrary D = 2n + 1 spacetime dimension. Our configurations are static and 
spherically symmetric, approaching at infinity a Minkowski spacetime background. Based on numerical re- 
sults for D = 5 and D = 7, 9, we conjecture that the presence of a CS term in the action allows for finite 
energy solutions for any D = In + 1, with n > 1. (The case D = 3 is special, since it requires the presence 
of a negative cosmological constant.) Most of the numerical results in this paper are for the D = 5, in which 
case, we provide a systematic discussion of the solutions in [28]. We have also discussed some results for 
D = 7, 9, which reveal some new features of the solutions. 

The paper is structured as follows: in Section 2 we present the general framework and analyse the field 
equations for an SO(D + 1) gauge group. In Section 3, the general features a consistent truncation of the 
general model for an SO(D— 1) x SO(2) gauge group are discussed. Numerical results for D — 5 and D > 5, 
respectively, are presented in Sections 4 and 5. We conclude with Section 6 where the significance of, and 
further consequences arising from, the solutions we have constructed are discussed. 

1 For the sake of completeness, one should mention that the picture is very different once one gives up the assumption of 
asymptotic flatness. For example, in anti-de Sitter (AdS) 3 + 1 dimensional spacetime, stable nA solutions have been shown 
to exist [151 ; there are also monopole and dyon solutions even in the absence of a Higgs field. As found in 1161 . some of the 
AdS nA solutions may provide a model of holographic superconductors. Also, the non-asymptotically flat nA solutions in |17| 
(with a dilaton field possesing a Liouville potential), have found interesting applications in providing gravity duals of J\[ = 1 
super- Yang-Mills theory. 
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2 The general model 



2.1 The action and field equations 

In odd spacetime dimensions, the usual gauge field action can be augmented by a (dynamical) Chern-Simons 
(CS) term. Restriction to odd dimensions follows from the fact that Chern-Pontryagin (CP) densities are 
defined only in even dimensions, and the CS density is defined formally in one dimension lower than the CP 
density. The resulting odd dimensional space is then interpreted as the spacetime on which the dynamical 
CS term appears in the Lagrangian. 

Such terms appear in various supersymmetric theories, the Af = 8, D = 5 gauged supergravity model |30j . 
[3"T] being perhaps the best known case, due to its role in the conjectured AdS / CFT correspondence. However, 
in this work we shall restrict ourselves to a simple EYMCS model, which does not seem to correspond to 
a consistent truncation of any gauged supergravity model. Also, our solutions approach asymptotically the 
Minkowski spacetime background. In the case of an Abelian gauge group in D — 5, a CS term leads to some 
new featureajj only for rotating black holes [33]. However, we shall see that for a nA gauge group, the CS 
term can affect the properties of solutions even in the static, spherically symmetric case. 

We consider the following action for the EYMCS model in D — 2n + 1 dimensions 



S = 



M \ lb7rG 



- £ym - K 



d D x C 



M 



CD) 
CS ' 



(2.1) 



where 



is the usual Yang-Mills lagrangian, with 



Fp U = d^A v - duAft 



■e[A^A v ] 



(2.2) 



(2.3) 



the gauge field strength tensor, while e and n are the gauge and the CS coupling constant, respectively. 
(The value of k is fixed in supersymmetric theories. However, in this work we shall treat k as a free input 
parameter. This has been also motivated by the study in [34] of the Einstein-Maxwell-CS system, which 
revealed a nontrivial dependence of the properties of the solutions on the value of k.) 

The definition of the Chern-Simons density on D— dimensional spacetime follows from that of the cor- 
responding Chern-Pontryagin density on D + 1 (even) dimensions. The latter is, by definition, a total 
divergence 

V • $7 = Tri F A F ■ ■ ■ A F 



n times 

in 2n— dimensions, and the CS density on a D = 2n — 1 dimensional spacetime is formally defined as one of 
the 2n components of the density ft. 

The CS densities Cqs thus defined are gauge variant. The explicit expressions of first three, in D = 3, 5 
and 7 dimensional spacetimes, are 



£g] = £ A ^TrMx 



C 



(5) _ Xpvpcr 



Tr < A 



^ 6 A^Ay 



FpvFpa — e F^ApAu + — e A^A v ApA a 



■•(7) Xfivparn 

-cs — £ 



Tri A 



p p p e p p A A p 2 F A F A 

pv x pa ± tk r pv ± pa^-r^K K pv- r *p ± ar-^n 



gfi 3 FpvApA a A T A K g^ e4 ApA l ,ApA a A T A H 



(2.4) 
(2.5) 

(2.6) 



2 Note that the situation can be different for charged magnetic branes, see e.g. the asymptotically AdSs Abelian solutions 
with a CS term in 1321 . 
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which are all manifestly gauge variant. Remarkably however, the Euler-Lagrange variations of these den- 
sities are actually gauge covariant. Indeed, these variational terms are expressed in gauge covariant form 
for arbitrary D = 2n + 1 as 



fi n _l_'\ \ ir A 1 lA l 2A I 3A 1 4---M2n-lM2n p TP TP 



M27i-lM2n 



Perhaps what is still more relevant in our case, where we restrict attention to static solutions only, is the 
fact that the CS densities (|2.4j) - (l2.6[) . etc., in that case reduce to a very useful form which can be expressed 
for the arbitrary D — 2n + 1 case. Working in a gauge such that dtA^ = 0, one can show that, up to a 
total divergence term (which we ignore here since we are only interested in the Euler-lagrange equations), 
the effective arbitarary n CS Lagrangian is 

4 2 s l+1) + e^*—- 1 *-Tr{ A F ilh F %3U . . . F t2n _ 1%2n \ . 



The field equations are obtained by varying the action (12.1[) with respect to the field variables g^ v ,A 



R„u - \g^R = 8irG T M!/! (2.7) 



9 

where 

(2.8) 



Tfiv = 2Tr|F MQ F^«? Q ' 3 - ~ 9fiv F a0 F a ^, 



is the energy momentum tensor. One can show that this tensor is covariantly conserved (i.e. VpT^ 1 " = 0) 
for solutions of the YMCS equations. 

In what follows, we shall be seeking to construct finite mass/energy solutions of the above equations. This 
is made possible by the fact that the energy density functional arising from the Lagrangian (|2.ip actually 
satisfies Derrick scaling by virtue of the presence of the CS term in it. This is because the CS term, which 
in D = 2n + 1 dimensional spacetime scales as L~( 2n+1 \ balances the Yang-Mills term which scales as L~ 4 . 
Specifically, in the 2n = D — 1 space dimensions, 2n + 1 > 2n > 4, for the cases of interest here, namely for 
n > 2. 

Inasfar as the CS term here plays the role of regularising the energy by providing the required Derrick 
balance, this makes it an alternative to employing YM higher order curvature terms |22[l25j for this purpose. 
The latter is of course more versatile since its use is not restricted, as in the CS cases, to 2n + 1 dimensional 
spacetimes. 

2.2 A spherically symmetric Ansatz 
2.2.1 The metric 

In D-dimensional spacetime, we restrict to static fields that are spherically symmetric in the D — 1 spacelike 
dimensions with a general metric Ansatz 

ds 2 = h(r)dr 2 + f 2 (r)dn 2 D _ 2 - f Q (r)dt 2 , (2.9) 

where r and t are the radial and time coordinates, while dQ 2 D _ 2 is the metric on the round (D — 2)-sphere 
(note that this Ansatz has still some freedom in the choice of the radial coordinate). 
The numerical work has been done for a metric gauge choice fzir) = r 2 and 

f 1 (r)=N(r) 1 f (r)=N(r)a 2 (r) where N(r)=l--^, (2.10) 



3 Strictly speaking Derrick scaling [24] applies only in flat space background. However, in practice it works also for gravitating 
configurations in asymptotically flat spaces 35] (at least in the spherically symmetric case). 
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the function m(r) being related to the local mass-energy density (as defined in the standard way) up to some 
Independent factor. 

Another convinient metric gauge choice used in the literature is 

/iW = -72- = Jff> Mr) = /(r), (2.11) 

corresponding to an isotropic coordinate system (the D = 5 exact solution discussed in Section 4.4 is found 
for this choice of coordinates). 



2.2.2 The YM fields 

The construction of a static static, spherically symmetric YM Ansatz in D — 5 spacetime dimensions leading 
to a nonzero CS term has been discussed in [36]. In what follows we present an extension of that result for 
a generic D = 2n + 1 case. 

There is some arbitrariness in the choice of the gauge group. The only restriction is that it should be large 
enough to accomodate for a static spherically symmetric Ansatz, with a nonvanishing electric potential- In 
D = 2n + 1 dimensions, the smallest simple gauge group supporting a nonvanishing CS term is SO(2n + 2). 
Here we shall take the SO(2n + 2) YM fields in one or other chiral representation of SO±(2n + 2). Our 
spherically symmetric Ansatz is expressed in terms of the representation matrices 



1 ± r 2r , 



\r a ,T p ] 



a, (3 = l,2,...,2n+2 



(2.12) 



T a = (Ti, Tm), with the index M = (2n + 1, 2n + 2), being the gamma matrices in 2n + 2 dimensions and 
T2n+3, the corresponding chiral matrijJl. We shall adopt a normalization convention such that 



2Tr<^ £ 



i(±) v (±) 



S aa 'S, 



0/3' ■ 



(2.13) 



Our construction of a spherically symmetric gauge field Ansatz is based on the formalism of A. Schwarz 
[3 7) . An alternative formalism, [38 , is familiar in the literature, but the calculus of |37j was found to be 
more convenient for the purposes of this work. 

Our spherically symmetric Ansatz for the YM connection — (Ao,Ai) is 



A = 



\{-{zx) M ^-X 



A = - 

e 



2n+3 v (±) 

Zj 2n+l,2n+2 (l 



6271+3 + 



lM 



-A 



2n+3 ~ 



(±) 



X i ^2n+l,2n+2 ( ) 



(2.14) 



(2.15) 



(with x 1 the usual Cartesian coordinates on R D ~ 1 and x° — t) in which the summed over indices M,N = 
2n + 1, 2n+ 2 run over two values such that we can label the functions (</> M , </> 2n+3 ) = tf>, (\ M , x 2n+3 ) = x 
and (A^ 1 ,A 2n+3 ) = A r like three isotriplets 0, x an( i A, all depending on the 2n dimensional spacelike 
radial variable r and time t. e is the two dimensional Levi-Civita symbol, while Xi — Xi/r (with XiXi — r 2 ). 

In what follows we are interested in configurations without a dependence on time. Then the parametri- 
sation used in the Ansatz (|2.14p - (|2.15[) results in a gauge covariant expression for the YM curvature 



4 For the most interesting case D = 5, this condition rules out the possibility of using the minimal non-Abelian gauge group 
SO (3) EJ. 

5 Thus, the fact that we are using an antihermitian representation for the SO(D + 1) algebra matrices leads to a factor of i 
in front of CqJ . 
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Ffxv — (Fij,F i0 ) 



Fi 



F 4 



in 



1 



M v(± 



£ r 
1 



2n+3 



:±) 



A/ 



\M1 



(±) 



DrX 



2n+3 



^1 S 



(±) 

2n+l,2n+2 



in which we have used the notation 

# r a = d r 4> a + s abc A b r 4> c , D rX a = d rX a +e abc A b rX c 



(2.16) 



(2.17) 



(2.18) 



), X = X" = (X M ,X 2 " +3 ), with 



as the 5*0(3) covariant derivatives of the two triplets (f> = <j) a 
respect to the 50(3) gauge connection A r = A^. 

After taking the traces over the spin matrices, it is convenient to relabel the triplets of radial function as 



— (aM 



), x=(x M ,X 3 ) and A r = {A^,A% 



— I am 



(2.19) 



with M = 1, 2 now. 

The triplet, A r (r), plays the role of a connection in the residual one dimensional system after the impo- 
sition of symmetry, and encodes the SO(D — 1) arbitrariness of this Ansatz. In one dimension there is no 
curvature hence it can be gauged away in practice |36) . 

However, finding solutions within the YM Ansatz (12.14[) . (|2. 15|) (which after setting A r = still features 
six independent functions), is technically a difficult task. A further consistent trucation of the general 



Ansatz is 



X 



0, leading to an EYMCS system with six unknown functions, four of them being 



gauge potentials parametrising the gauge field, and, two metric functions. Indeed, the two gauge functions 
suppressed are redundent and would only be excited in an eventual stability analysis of our solutions. To 
make connection with previous results on EYM solutions, we shall note 



(p 1 (r)=w(r), 3 (r) = w(r), x 1 (r) = V(r), x 3 (r) = V(r), with A r = 0. 



(2.20) 



The Ref. [36] gave numerical evidence for the existence of solutions within the above Ansatz for the case of 
D = 5 (i.e. an 50(6) gauge group) and AdS asymptotics. Some of the features discussed there are generic. 
For example, the resulting system has some residual symmetry under a rotation of the 'doublets' w(r), w(r) 
and V(r), V(r) with the same constant angle u (e.g. w — > w cos u + w sin u etc.). Note that for configurations 
with w(r) — V(r) — the gauge potentials are unvariant under the "chiral" transformations generated 
by X^ri+i 2n+2' ^ ne configurations with w(r) — V(r) = instead change just by a sign under the same 
transformations. Also, this Ansatz is invariant under the parity transformation 4> a —¥ —cf> a , x" ~* ~X a - 
A further simplification of the YM Ansatz consists in taking 



w(r) = V(r) = 0, 

which is a consistent truncation, SO(D — 1) x 50(2), of the general Ansatz. 



(2.21) 
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2.3 The equations 

The truncated YM Ansatz (|2.20[) together with the generic metric element (|2.9p leads to the following set 
of EYMCS equations: 

fK-&-$T + ^(/o/2 + 2/oA - ^f) -^r((D- 2)(2D - 7)h(Vw - Vwf (2.22) 
2/0 2/1 2/2 2/2 4/2 V 

+2(2D - b)h{V' 2 + V' 2 ) + 3(-P- 2 )(^-3)/o/i (w2 + ^ 2 _ 1)2 + 2(jD _ 2)Uw ,2 + = Q) 

72 / 

+ (J - 3)/ 2 2 _ _ 3)/i + a2 f h {v a + y l2) _ {D _ 2){w a + ~, 2) (2 23) 
Jo 2/2 V 70 

H TTT ( w (Vw-Vu;) =0, 

2/2 2/0 / 

- H + - - 3)/i + y (P - 2)(^' 2 + w' 2 ) + ! f{V' 2 + V' 2 ) (2.24) 

( D { v w - vtif + ( ^- 2) if- 3)/l (^ + J- I) 2 ) = 0, 



2/0 ^ ' 2/; 

// , l//o / 1 1 (-^ — 4)/2- , /l fV/f> t^~\ (-C — 3)/i . 2 2 . . . 

w +^(-r~^-H 7 )w + — ^(Fw - Fuj) w(w +w -1) (2.25) 

^ jo h h 2jo J2 

m ^-^' (^ + - 2 - D ( "- 5)/2 (V + - 2 - 1 ) v " + P - 3 )(^- - ^)-') = 0. 

W ~ 2j/ 2 V 71 70 v / 

-// . l//o / 1 1 (-^ ~ 4)/2, , /l T/ / T/ ~ fV \ {D — 3)fl 2 1 ~2 i\ ; nc \ 

w +r(-T---T-H 7 )iu + — K (Fiy - Viu) w(u> +w -1) (2.26) 

^ 70 7i 72 2/o 72 

' m T^Wrr-r ^ + w2 - 1)( " 5)/2 ( ^ + " 2 " ^ + ( ^ ~ 3)( ^ ~ VW)W ') = °' 

(,^7 _ ^J/2 V 7170 \ / 



p+ ( 2 ir fi 44) f, - 2 if ! i f '- ff i (2 - 28) 

\ 2/2 2/0 2/1 y 2/2 
-«(£> 2 - 1) sfhhfi- DI \v? + w 2 1) (D ~ 3)/ V = 0, 

where a prime denotes a derivative with respect to r. Also, to simplify the expression of the above relations, 
we note 

o 16irG .„ „. 

° - (7^2)^ (2 - 29) 

and absorb a factor of l/e^ D ~ 3 ^ 2 in the expression of k. 

After fixing a metric gauge, the eq. (|2.23[) becomes an Hamiltonian constraint. There is also a constraint 
equation for the gauge fields, 

D-2 

f 2 ~ Q-4 / t n 

72 /TAX/' T/I^-l 1 ^7-1 r,\fT- /70/ - / 



(W - VV) + (D- 2)/ 2 2 \ ^{ww' ~ ww') - k(D z - l)(Vw - Vw){w z + w z - 1)~ = , (2.30) 

v/0/1 y 71 

which originates from the variational equation for A r (one can show that (|2.30p is a first integral of the eqs. 
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Also, in what follows, we shall restrict to a dimension D > 5. The case D — 3 should be discussed 
separately, since the existence of physically interesting solutions requires the presence of a cosmological 
constant. 

2.4 The asymptotics and a truncation 

Numerical evidence for the existence of asymptotically AdSs solutions within the full Ansatz (|2.20p was 
given in Ref. [3B]. However, it seems that the presence in that case of a negative cosmological constant was 
crucial in arriving at that resulfl In the asymptotically fiat case, we could not find such solutions (with 
four essential functions) but only configurations within the restricted SO(D — 1) x 50(2) Ansatz (I2.21|) . 

Although we do not have a rigurous proof of that, some analytical indications in this direction come 
from the study of the asymptotics of the general solutions close to the horizon and at infinity. Here it is 
convenient to use the metric Ansatz (|2.10[) . with two functions N(r) and a(r). For black hole solutions, the 
horizon is located at r = > 0, with N(rh) = and a(rh) > 0, while A'(r^) > in the nonextremal case. 

Then we suppose that all functions admit the following behaviour asr-> r/>: 

oo oo 

w(r) = cosC/i E w k( r - r h) k , w(r) = sinf/i E*fc( r - r h ) k , 

k=0 k=0 



V(r)=cosU 1 J2v k {r~r h ) k , V(r) = sin U 1 ]T i k (r - r h ) k , (2.31) 

fe=0 fc=0 

CO oo 

N(r) = Y / N k (r- r h ) k , a(r) = E a k (r - r h ) k , 



k=l k=0 

which satisfy the regularity condition wV — wV — > as r — > r^. The cofficients w kl w k , V k , m k and a k 
are computed order by order by substituting this expansion in the field equations. It turns out that the 
only free parameters are cr(r^), w(r^) and v\. Moreover, we have verified that, at least up to order fouiQ 
w k /w k = V k /V k = 1. 

Interestingly, a similar analysis for large values of r leads to the same conclusions. Here we suppose the 
solutions admit a power series expansion with 

w(r) = cosU^ — jp w(r) = sinU^—^-, V{r) = cosZJ^ V(r) = cost/^ -|, 

k=0 fc=0 fe=0 k=0 

oo . , oo 

k=D-Z k=l 

After plugging this expansion in the field equations, we have found that, W k = W k and V k = V k , at least up 
to order D + 5. Moreover, the only free parameters in the above expressions are W2 and v%. 

This result, together with the corresponding one for the near horizon expansion (|2.3ip strongly suggests 
that the functions iv(r), w(r) and V(r), V(r) have a constant ratio for any r > for any physical solution of 
(|2.22p - (|2.28p . Although we do not have a rigurous proof, this conjecture has been confirmed by our numerics 
and all black hole solutions we have found have in fact only two essential gauge functions^. This applies 
also for asymptotically flat solitons {i.e. without an event horizon), in which case we have also failed to find 
asymptotically solutions within the general Ansatz (|2 . 20|) . In asymptotically AdSs spacetime, the asymmetry 
between w, w and V, V explictly appears in the large-r behaviour, being introduced by the cosmological term, 
see the results in Section 2 of Ref. [36] . 



6 This is not an entirely surprising result. We recall that already in D = 4 dimensions and a gauge group SO (3), the presence 
of a negative cosmological constant A leads to some new qualitative features |15l . In particular, A < allows for EYM static 
solutions with a non- vanishing electric potential, which is not the case for asymptotically flat configurations |39| . 

7 Beyond this order, the involved relations were too complicated to deal with. 

8 Although we could construct D = 5 black hole solutions within the general ansatz J2 . 20ft , it turns out that, within the 
numerical accuracy, the ratios w(r) /w(r) and V(r)/V (r) were in fact always constant. 
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Then, for the remaining of this work we shall deal with the case of solutions within the restricted 
SO(D — 1) x 5*0(2) Ansatz with two esential functions: a magnetic potential, w(r), and an electric one, 
V(r). 

3 The SO(D - 1) X S0{2) model. General properties 
3.1 The equations and scaling properties 

The equations of the model simplify drastically for the truncation (|2.2ip . Their form within the metric 
parametrization (|2.10p . reads 

W +[ - + N + a )W (D-2) aNrD-i V + r *N ~ 



r 



(3.1) 



V ' 2 ■ ^- 2 ^' 2 + (D-2)(D-3)( 1 -^ 



2 V cr 2 r 2 y 2r 

a' = a 2 p-2)^-, 
2r 

the gauge constraint (|2.30p vanishing identically. These equations can also be derived from the effective 
action 



S e ff — I dtdr^ am! — —a 2 



,. X AW 2 V' 2 (D-2)(D-3) , 9x9 . 
' -[( D ~ 2 )—2 - + - £i Uil-w 2 ) 2 ) (3:1) 

- 2k(D 2 - l)V(w 2 - l)^w' 



(note that, as required, there is no coupling with the geometry for the term proportional with k). 

A generic feature of the YMCS model within the SO(D — 1) x SO (2) truncation is the existence of a 
first integral for the electric potential V(r), 

V = (4 + (D 2 - 1)kF(w)) , (3.3) 



where P is an integration constant (we shall see that, for globally regular solutions, it's value is fixed by k). 
The function F(w) has the following general expression in terms of the hypergeometric function ^Fy. 

FW = (-1)» (W) 2F 1 {\^,\^)w, (3.4) 
its explicit form for several dimensions being 

F(w) = w, for D = 3; F(w) = —w + — w 2 , for D = 5, 

o 

2 w 5 

= w w 3 H , for D = 7. 

3 5 

3 1 

f(w)=-w + w 3 w 5 + -w 7 , for D = 9. 

5 7 

One should also note that the eqs. (|3.1[) together with the first integral (|3.3p are invariant under the scaling 



and 



r 



Ar, m ->• A- 3 m, er ^ er, w -> w, -> V/A, P -> X D 2 P, and a -> Aa, k ^ X D 4 k, (3.5) 
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with A an arbitrary positive parameter. There is also a second scaling symmetry of the equations (|3.1[) 

V -> XV, o — y \a, (3.6) 

together with t — t/A, all other variables remaining unchanged. This symmetry is lost after setting er(oo) = 1 
as a boundary condition. 

The last symmetry of the equations of the model consists in simultaneously changing the sign of the CS 
coupling constant together with the electric or magnetic potential 

k —t —k, V — >• — V, or k — y —k, w —> —w, (3-7) 

(the first integral (13.31) implies also P —> —P in the first case). In what follows, we shall use this symmetry 
to study solutions with a positive k only. 



3.2 The behaviour at infinity 

Unfortunately, it is not possible to find an exact solution of the equations (|3.ip with a nontrivial magnetic 
gauge potential w(r), except for a special value of k in D = 5 dimensions. However, one can write an 
approximate form of the solutions as a power series with a finite number of undetermined constants, both 
at infinity and at the horizon/origin of the coordinate system. This analysis allows us to obtain some 
information on the possible global behaviour of solutions. 

In deriving the asymptotic form of solutions as r — > oo, wc assume that the spacetime approaches the 
Minkowski background at infinity, while the configurations have a finite ADM mass. Then we arrive at the 
following expansion of the solutions at r — > oo: 

1 O 2 a 2 (D-3) 2 J 2 

m(r) = M --a\D-3)^ + ..., a(r) = 1 - +..., (3.8) 

w (r)=±l--^ + ...,V(r) = V --^ + ... . 

In the above relations, J, Mq, Vq are parameters given by numerics which fix all higher order terms, while 
Q is a constant fixing the electric charge of the solutions, 

' % + (D 2 - 1)kF(±1)) , (3.9) 



" D-3 
where 



F{-1) = -F(l) = (-1) — 



r(V) 

2T(f ) ' 



The set (I3.8P of boundary conditions is shared by both globally regular and black hole solutions. 



3.3 Soliton solutions: the expansion at r = 

These are perhaps the simplest possible solutions of the system (|3.1[) and can be viewed as higher dimensional 
generalizations of the Bartnik-Mckinnon solutions [13], though dressed with an electric charge. The most 
striking feature here is that the electric charge of the solitons is fixed by the value of the CS coupling constant. 
Technically, this results from the fact that the term P/a 2 + (D 2 — 1)kF(w) in the first integral (|3.3p should 
vanish as r — > 0. Since w(0) — 1 for regular solutions, the parameter P is fixed to be 

P = 2 a 2 (-1)^ K ^^^. (3.10) 
One finds e.g., P — 16a 2 n, — 128q; 2 k/5 and 256a 2 n/7 for D = 5, 7 and 9, respectively. 
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Then, for globally regular soliton-type solutions, the electric charge parameter which enters the far field 
expression (|3.8[) is fixed by 

c - cM _ j; 4(-D^yyr(^) 
Q-Q -« (Z) _ 3) r( ^ } ■ 

Its expression for the values of D considered in numerics is = 16k, -64k/5 and 256k/21 for D = 5, 7 
and 9, respectively. 

Also, one finds that the globally regular configurations have the following expansions near the origin 
r = 0: 

w (r) = l - br 2 + 0(Y 4 ), V(r) = (-l) i ^2^6- 2 ^ i (£) + l) K a r 2 + 0(r 4 ), (3.12) 
m(r) = (£> - 2)a 2 & 2 r D ~ 1 + 0(r D ), a(r) — a + (D ~ 2)a 2 b 2 a r 2 + 0(r 4 ). 

The free parameters here are b — —^w"(0) and (Jq = cr(0). The coefficients of all higher order terms in the 
r — > expansion are fixed by these parameters. 



3.4 Black holes: the near horizon solution 

For the solutions in this work, the event horizon is located at a constant value of the radial coordinate r = r/j , 
with N(rn) — 0. In this case, one can write also an approximate form of the solutions near the horizon, as 
a power series in r-r^. In the nonextremal case, the first terms in this expansion are 

m(r) =r h + mx{r - r h j + . . . , u{r) = a h + — (r -r h ) + 

w(r) = Wh + wi(r - rn) + . . . , V(r) =vi(r- r h ) + (3.13) 



where 



= (P + *\D 2 - 1)«FK)) , mi = ir°~ 2 ( * + (D - 2){D - 3) (1 f 



(£>-3KK-i) (tf-iy^K-iiT ^ 2 P-2) 2 

Wl= r h (D-3-^) + P-2K(- TOl + (^-3).r 4 )' ^ = ^— ^ 

h 

while AT(r) = N'(r h )(r - r h ) + . . . , with N'(r h ) = (D -3-m 1 /r h J ~ 4 )/r h . The obvious condition N'(r h ) > 
implies the existence of a lower bound on the event horizon radius r m ; n , for given values of Q and n. The 
only free parameters in the expansion above are ah and Wh- In a numerical approach, their values are found 
by matching the near horizon form of the solutions (|3.13l) with the asymptotic expansion (13.8[) . 

As the function N(r) develops a double zero at the horizon, i.e. N(r) = N 2 (r - rh) 2 + . . . , 

and the black holes become extremal. We shall see that such solutions exist indeed, emerging as limiting 
configurations of a branch of nonextremal black holes. The near horizon expansion is more constrained in 
this case. Supposing w\ ^ 1, one finds that the event horizon radius rh is the largest positive solution of the 
equation 

the value of the electric charge parameter being also fixed (with w(oo) = ±1): 

q - r ff-y^'g ^*.) - n„» - a . (3.i6) 



(D 2 -1)tt* \{D-2)(D-3) 



w h r i 
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Also, the coefficient N 2 in the leading order expansion of the metric function N(r) is 



No = 



(£>-3)(D-4) 



o r 2 



a 2 {D - 2)(D - 6)(w 2 - l) 2 / 2{D - 2) 2 {D - i)r h 



2(13-4) 9 



4(L> - 4)r 2 



(D-6)(L> 2 -1)k 2 (w 2 -1) 



D-3 



- 1 



(3.17) 



The parameters w\ , ui and a\ have a similar expression as those found in the nonextremal case (there one 
should replace the expressions (|3.15|) and (|3.16j) for r/j and Q, respectively). 

For D = 5 and Q > Q^, we have found numerical evidence for the existence of a different type of 
extremal black holes, with w(rh) = 1. As r — > r/,, these solutions have basically the same leading order 
expression as an extremal RN black hole. However, the next to leading order terms in the near horizon 
expansion exhibits non- integer powers of r — r/,. . These special D = 5 configurations will be discussed in the 
Section 4.3 below. 



3.5 An AdS 2 x S D ~ 2 solution 

As expected, the near horizon structure of the extremal solutions with w 2 7^ 1 can be extended to a full 
AdS2 x S D ~ 2 solution of the field equations. This is a new exact, essentially nA solution to the EYMCS 
field equations, with 

dr 2 

ds 2 = Vl (— -r 2 dt 2 ) + v 2 cM 2 D _ 2 , w(r)=w , V(r)=qr, (3.18) 



where 



(D-2)0D-3)t>n4 D ~ 6)/2 w 



H (D 2 -l)n (W 2 - l)(£>-5)/2 

The AdS radius satisfies the relation 



(3.19) 



a/Vn o^Afr, ^ , am o\ / 2 ^ c n , 2(J> - 2) 2 (D - 3)w%v?- 4 

(D 2 - 1) 2 k 2 (w 2 - 1) D " 3 



vi = 8v 2 2 (D - 3)(4(U - 4)«a + c^(£> - 2)(w 2 - 1) 2 (6 — D + v / v 2 ' ' 2 ) , (3.20) 



where the size of the S D 2 part of the metric results as a solution of the equation 

1 2/ n 9 v 2 n a A , 2(£>-2)(£>-3) iu 2 A 
W2 = -a (D - 2)K - 1) ^1 + p 2 _ 1)K 2 (w g-ip-3 ^ J - (3-21) 

being a function of wq only. Unfortunately, one can write a simple solution of the above equation, except 
for D = 5. In this case, the general solution reads 

1536a 2 K 4 (wg - l) 2 48a 2 n 2 (w 2 - l) 2 _ 32V3aK 2 w (w 2 - 1) 

Vl ~ (64k 2 - a 2 w 2 )(M K 2 + a 2 w 2 )' ^ " 64 K 2 - a 2 w 2 ' Q ~ ^ 6 4« 2 - ah^^K 2 + a 2 W 2 ) ' ( ' } 

One can see that, for any D, the properties of the AdS 2 X S D ~ 2 solution are uniquelly specified by the 
constant wq. Note that for Wq 7^ 1, the nA magnetic gauge field is nonvanishing; the case wq = ±1 is special 
and describe the near horizon geometry of the extremal RN solutions. 

Although finding local solutions in the vicinity of the horizon does not guarantee the existence of the global 
solutions, the above result provides an argument that the EYMCS system is likely to present asymptotically 
fiat, extremal black hole solutions. In D = 5 and D = 7, 9, this is confirmed by our numerical and analytic 
results. If there are extremal black holes in the bulk, the parameter q in (|3 . 1 8|) is related to the bulk charge 
parameter Q via 

vi 



(D-2)/2 

v 2 



((£> - 3)Q - (D 2 - 1)k(F(w ) + F(±l)) , (3.23) 



with to (00) = ±1, the allowed values at infinity of the bulk magnetic gauge potential. 

It would be interesting to consider these solutions in the context of the attractor mechanism and to 
compute their entropy functions (for a discussion of Sen's entropy function for D = 5 supergravity models 
containing Abelian Chern-Simons terms, see e.g. |40) ) . 
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3.6 Relevant parameters and global charges 

The only global charges associated with the solutions are the mass M. and the electric charge Q 

M= ( d -l V ^ M ° , Q= <P-Wn-*Q (324) 
107rG g 

where Vb-2 = 2tt^ d ^ 1 ^ 2 /T((D — l)/2) is the area of the unit D — 2 sphere. The mass M. is the charge 
associated with the Killing vector d/dt; the electric charge Q is associated with the U(l) gauge symmetry 
generated by S^n+i 2n+2- For the boundary conditions in this work, the electrostatic potential difference $ 
between the horizon and infinity is fixed by the value at infinity of the electric potential V(r), as read from 
(|3.8[) . In a thermodynamical description of the system, $ corresponds to the chemical potential, 

$ = — . (3.25) 

e 

The Hawking temperature and the entropy of the black holes are given by 

T H = ^-a{r h )N'{r h ), S = with A H = U D _ 2 rf ~ 2 . (3.26) 

47T 4(jr 

The solutions should also satisfy the First Law of thermodynamics, 

dM = T H dS + <5>dQ. (3.27) 

In the discussion of the black hole thermodynamics, we shall restrict to configurations in a canonical ensemble, 
the relevant potential being the Helmholz free energy 

F[T H , Q}=M- T H S. (3.28) 

As usual, in practice, it is convenient to work with quantities which are invariant under the rescaling (|3.5p 
(note that in the numerics we set a = 1). For solutions in a canonical ensemble, we normalize the global 
quantities with respect the charge parameter Q and define e.g., the dimensionless quantities 

/ = 757-2' tH = T H Q^ D - 3 \ a H = ^ L r and j = L. (3.29) 
Q/9 Qt^ Q 

One should note that there is no conserved quantity associated with the parameter J which appears in the 
large-r asymptotics Q3.8p . Also, J does not enter the First Law p.27[) . 

In addition to the mass and the electric charge, there is another global charge which has a topological 
origin. This is the volume integral V of the topological density calculated from the "magnetic" components, 
Fij, namely the Chern-Pontryagin (CP) density defined on the (D — 1)— dimensional space dimensions with 
Euclidean signature. The correct expression for this CP density must take account of the fact that the 
gauge group for is SO{D — 1), and not H one or other of the two chiral algebras SO±(D — 1). Using the 
spherically symmetric components of the curvature Fij given by (|2.16|) . this charge density is calculated to 
be 

e lllwll ... lD ^ lD _ x Tr {^ D>D+1 F^ . . . F iD ~^ } ~ --^ (w 2 - lJ^rV . (3.30) 

Usually, for spherically symmetric soliton solutions, the integral of the above quantity is suitably normalised 
such that it yields unit "magnetic" CP charge, V = 1. However, noticing that (w 2 — l)^~w' is just the 
derivative of the function F(w) which enters the first integral (|3.3|) . it is more interesting to use a nonstandard 
normalization and to define the "magnetic" CP charge V directly as the integral of (|3.30l) . Then, taking into 
account the boundary conditiont0 (13. 8[) . (|3.12[) . it follows that the soliton solutions exhibit an interesting 
connection between the electric charge and the "magnetic" CP charge 

Q = kV. (3.31) 



9 Note that the SO(D-l) curvature consists of two SO±(D — 1) curvatures, each contributing ±1 CP charge in the spherically 
symmetric case. If the appropriate factor of = Id m 113.301 1 is not accounted for in the CP density, these two charges 

will cancel each other out. This is discussed in detail in the context of the SO(D) monopole on TR D in [4T] . 
10 Note that we take ui(oo) = —1, which was the case for all solutions we could find numerically. 
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3.7 The Reissner-Nordstrom solution 

The Reissner-Nordstrom (RN) solution plays an important role in what follows. Thus, for completness, 
we briefly discuss here its basic thermodynamical properties. This solution is recovered for w(r) = ±1, in 
which case the magnetic components of the field strength vanish identically. Then the configuration becomes 
essentially Abelian and one finds the following exact solution 

m(r) = Mo - a2( ^~_3 )Q2 , cr(r) = l, u>(r) = ±1, V(r) = V - (3.32) 

with Vq a constant which is usually chosen such that the electric potential vanishes on the horizon, i.e.,Vo = 
Q/r®~ 3 . The mass M., electric charge Q and chemical potential $ follow automatically from (|3.24l) . f!3.25[) . 
The Schwarzschild-Tangerlini vacuum black hole corresponds to the case Q = 0. 
The RN solution has an outer event horizon at r — r^, with 

1 , \ l/P-3) 



r h = [^(M + ^Af 2 - 2a 2 (D - 3)Q 2 )J , (3.33) 

whose existence imposes an upper bound for the electric charge for a given mass. The Hawking temperature 
and the entropy of this solution can be written in terms of r^, Q (which are the parameters used in numerics) 

as 

With these definitions, one can easily verify that the First Law (|3.27[) is indeed fullfilled. In addition, the 
RN black holes satisfy the Smarr law 

M = ^^T H S + <S>Q. (3.35) 

The relation (|3.34p shows that, for a given Q, there is a minimal value of the event horizon radius at 
which an extremal black hole is approached, 



</• >r { r n) =(\{D-i)a 2 Q 2 y D ^ . (3.36) 

Thus the Hawking temperature vanishes in this limit, while the event horizon area approaches a minimal 
value, with 

a' mm) = i(D-3) Q 2 ^, (3.37) 

After expressing as a function of Ah according to f|3 . 26[) . one gets the following relation between the 
reduced variables t# and an- 

Unfortunately, one cannot invert this relation to get S(Th, Q)- However, one can see that, for a given Q the 
solutions exist only for < Tjj < j , ^" csx ) ) the maximal value of the Hawking temperature being 

T { ™ ax) = r }_ 7) (£>-3)^ J 5T(2£)-5)~^^ L_ ; (3.39) 

2(20-6)71- (aQ)°-3 
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the entropy at this turning point being 



5D-14 
2 2(D-2) 



P — 2 
D — 3 



G \ a v /(L>-3)(2D-5)Q 



(3.40) 



It is also of interest to express the free energy of the RN solution as a function of Th, Q- The only relation 
we could find in this case reads 



/ = 
I-/ 2 



1 



(D-3) 



D-3 



Vd-2 ( , L 2(X>-3)(2£>-5)V^_ 2 1 ^ 1 



ig- 3 (£> - 2)2 2C - 7 7r c - 3 a 2 
a 2 (,D-2) 2 



1 + 1 1- 



a 2 (L>-2) 2 



(3.41) 



2(2i?-5) 2 ( J D-3)^_ 2 



1-4/1 



2(£>-3)(2£>-5)yg_ a 1 



a 2 (L>-2) 



,f 2 



D-3 



A study of the above relations shows the existence of two branches of Abelian solutions, each with different 
thermal properties. (The generic picture here is dimension independent; then the well-known D — 4 result 
in [42] applies also in higher dimensions.) For a fixed Q, there is first a branch of large black holes whose 
entropy decreases with Th, which therefore are unstable. This branch stops in a critical configuration with 
a maximal value of Th given by ()3.39j) . where a secondary branch of small black hole emerges. This branch 
has a positive specific heat and ends in an extremal configuration with an event horizon area given by (|3.37j) . 
These features are exhibited in Figures 7, 8 (see the RN curves there). The picture for D > 5 is qualitatively 
the same. 



3.8 The issue of perturbative static solutions around Reissner-Nordstrom black 
holes 

Since the RN black hole is a solution of the model for any D, one might expect the existence of a branch of 
nA solutions connected with it. Such solutions, if they exist, would emerge as static perturbations around 
the RN background. However, as we shall argue, this is the case for D = 5 only. 

Suppose that there is a perturbative solution of the equations (|3.1|) around the RN black hole, 

m(r) = m (r) + emi(r) + . . . , <r(r) = 1 + e<ri(r) + . . . , (3.42) 
w(r) = ±1 + eWx(r) + . . . , V(r) = V (r) + eVi(r) + . . . , 

with e a small parameter. In the above relations, mo(r), Vo(r) are the functions which enter the RN solution 
(13321) . 

After substituting (|3.42p in the equations (|3. 1[) . one finds that to lowest order, the equation for W\{r) 
decouples. This equation is the only relevant one, it's general-D expression being 

(r D - 4 NW[Y = 2((D- 3)r d - 6 ± 8nV^5 D ^) W u (3.43) 

(with N = 1 — mo(r)/r D ~ 3 ). Then it turns out that the case D = 5 is special, since the CS term gives 
a nonzero contribution to the W^-equation only for this dimension. The perturbation W\{r) starts from 
some (arbitrary) nonzero value at the horizon and vanishes at infinity, in order to be consistent with the 
asymptotic behaviour (|3.13[) . (|3.8[) . However, one can show that for D > 5, there is no solution of (|4.14l) 
that satisfies this asymptotic behaviour. To prove that, we rewrite the eq. (|3.43|) in the equivalent form 

~(r D - A N(W?yy = ^r D ~ 4 NW[ 2 + 2{D - 3)r°- 6 W? , (3.44) 

recalling that D > 5 now. Then, after integrating from to infinity, one finds that W\ necessarily vanishes 
identically^. The same argument applies when considering higher order terms in the expansion (|3.42p . 

^Basically, the r.h.s. of 113.4411 is greater or equal to zero, while N{r h ) = and Wi(oo) = 0. 
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Therefore we conclude that, for D > 5, the RN solution is stable with respect to nA perturbations within 
the considered EYMCS model 0. 

For D = 5, a similar reasoning implies that the "mass" term 1± in the Eqn. (|3.43p should necessarily 
be negative in the vicinity of the horizon. Then, one finds kQ > for w(po) — — 1 and kQ < for w(oo) = 1. 



4 The results in D = 5 
4.1 Numerical methods 

The scaling transformation Q3.5P can be used to fix an arbitrary valu d 1 " 3 ] for a. The usual choice is a = 1, which 
is what we employ for all solutions in this work. This fixes the EYM length scale L = y/l6nG/(Jp — 2)e 2 ), 
while the mass scale is fixed by fj, = L^ D ~ 3 ^ 2 /G. All other quantities get multiplied with suitable factors of 
L. 

To control the quality of the numerical results, we have performed some of the calculations with two dif- 
ferent methods, finding excellent agreement. First, the equations (|3.1[) were solved with suitable boundary 
conditions which result from (|2.31|l . (12.32ft using a standard solver [43] ■ This solver involves a Newton- 
Raphson method for boundary-value ordinary differential equations, equipped with an adaptive mesh se- 
lection procedure. Typical mesh sizes include 10 3 — 10 4 points. The solutions in this work have a typical 
relative accuracy of 10~ 7 . In this approach, the value of the electric potential at infinity Vq is fixed, the 
electric charge resulting from numerics, i.e., the configurations are in a grand canonical ensemble. (The first 
integral (|3.3|) has been used to verify the accuracy of the solutions.) 

In addition to employing this algorithm, families of solutions with a fixed electric charge were constructed 
by using a standard Runge-Kutta ordinary differential equation solver. In this approach we evaluate the 
initial conditions at r = r% + 10~ 5 , for global tolerance 10~ 12 , adjusting for shooting parameters and inte- 
grating towards r — > oo. In this case the electric charge Q was fixed via the equation (|3.3p . the electrostatic 
potential Vo resulting from numerics. We have confirmed that there is good agreement between the results 
obtained with these two different methods. 

Also, for both approaches, we have restricted our integration to the (physically more relevant) region 
outside of the horizon, r > r^. 



4.2 Perturbative solutions: an instability of the RN 5 black hole 

The branch of D = 5 nA solutions emerges as a perturbation of the RN black hole. (In what follows, we 
shall suppose without any loss of generality that the RN black hole has a positive electric charge, Q > 0.) 
The perturbative solutions are found by solving the Eqn. (|3.43l) . which for D = 5 reads 

r(rNW[)' -A(l±^-)Wx = ti , (4.1) 

Q 2 + r 4 O 2 

where N — 1 — 2 2 h + . Although this linear equation does not appear to be solvable in terms of known 

h 

functions, one can construct an approximate solution near the horizon and at infinity. For a vacuum choice 
w = — 1, one finds that, as r — > 

W l{ r) = W h + 2Whr t 8 T rl) ^-^ + --- ^ 

with all higher order coefficients fixed by Wh- Because (|4.2p is linear, one can take Wn = 1, without any loss 
of generality. 



12 At the cost of replacing the usual Yang-Mills term F 2 by F 2p |51| . one can expect that a brach of the static nA solution 
emerges as a perturbation around the corresponding RN-type background, in the D = 4p + 1 case, also. 

13 In principle, one can use 113.51 1 to fix instead the value of the CS coupling constant k. However, this choice is less interesting. 
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Figure 1. The profiles of typical solutions of the D = 5 perturbation equation (|4.1[) are presented as a function of 
the radial coordinate r for two different values of k. The values of the charge parameter Q for the corresponding 
critical Reissner-Nordstrom solutions with Th = 1 are also shown. 




0.5 1 1.5 2 2.5 3 

K 



Figure 2. The shape function U(k) = Qjr\ which gives the unstable Reissner-Nordstrom solution where a branch 
of non-Abelian configurations emerges is shown as a function of the CS coupling constant k. The inlet shows the 
scaled free energy and Hawking temperature of critical Abelian solutions. 



At infinity, the only reasonable asymptotics reads 

J ( 2 J(Q 2 - AnQr 2 



W 1 (r) = ^+ ^ + -> (4 - 3) 



in terms of a free parameter J. The solutions interpolating between (|4.2p and (|4.3[) are constructed numeri- 
cally, typical results being shown in Figure 1. 

It turns out that such perturbative solutions do not exist for arbitray values of (77,, k). Restricting to 
solutions of (|4.1[) with monotonic behaviouJ^l of Wx(r), we find that for given (r^, k) for which such solutions 
exist, these pertain to a fixed value of Q. This value results from the numerics. 

The existence of these configurations can be understood as follows. For w(oc) = —1, the second term in 
(|3.43[) shows the existence of an effective mass term n 2 for W\ near the horizon, with ~ 1 — SnQ/r 2 ; all 



*Note that there are also solutions of 114.11 1 where Wi(r) has nodes. 
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solutions we could find have jj? < 0. Then it is also convenient to introduce the dimensionless function 



U(k) = ^f, 
rt 



(4.4) 



which uniquelly fixes the parameters of the critical RN solution. One finds e.g. 



VU(1 - U 2 




(4.5) 



The shape of the function U (k) for the fundamental set of solutions of (|4.1[) . for which W\ (r) are nodeless, 
is shown in Figure 2. The inlet there shows the scaled free energy and temperature of the critical RN solutions 
where static linear nA perturbations arise, as a function of k. 

The function U(k) depends monotonically on the CS coupling constant k. As k — > 1/8, one finds U = 1 
i.e., the corresponding RN solution becomes extremal, while the D = 5 RN black hole with the maximal 
value of the temperature (|3. 39[) is unstable for n ~ 0.55. Also, one finds that the function U decreases along 
the branch of large black holes, with U ~ l/4re for large k. 

No physically reasonable solutions of eq. (|3.43l) are found for k < 1/8, or for perturbations of the form 
w(r) = +1 + eW(r), in which cases the effective mass for W is always real (we recall we take Q > 0, k > 0). 

4.3 Nonperturbative black hole solutions 
4.3.1 General properties 

The instability discussed above signals the presence of a symmetry breaking branch of nA solutions bifur- 
cating from the RN black hole. 

Indeed, our numerical results provide evidence for the existence of finite mass black hole solutions of 
the EYMCS system with nontrivial magnetic gauge fields outside the horizon. These solutions smoothly 
interpolate between the asymptotics (|3. 13|) . p.8[) (all D = 5 configurations have u>(oo) = — 1). 

On the basis of analytical and numerical results we have a pretty clear picture of the behaviour of the 
D = 5 EYMCS black holes, this case being studied in a systematic way. The properties of the solutions 
depend on the value of the CS coupling constant k. For a fixed r^, the value of the electric charge parameter 
Q is also important, some basic features of the solutions depending on whether Q is less or greater than 



As a general feature, we could not find configurations with multinodes of the function w{r). Heuristically, 
this can be understood as follows: the existence of such configurations would imply w(r) = as the limit of 
multinodes. However, such solutions would have infinite energy, with m(r) = 3/21ogr + Mo, <x(r) = 1 and 
V(r) = 0, which is not compatible with the boundary conditions in the present work. 

More importantly, since our solutions emerge as perturbations of RN black holes, we notice the existence 
of configurations without nodes of the magnetic potential w(r). We shall see that some of these solutions 
are stable. Also, different from the case of other asymptotically flat hairy black holes with nA fields |12j . 
w(r) may take values outside the interval [—1, +1]. 

Moreover, it is possible to find more than one solution for the same value of (k, Q, r^). In this case, apart 
from configurations with a monotonic behaviour of w(r), there are solutions solutions with local extrema 
of the magnetic gauge potential, see Figure 3 for such an example. However, it is likely that the solutions 
with local extremal are always thermodynamically disfavoured because spatial oscillations in w increase the 
total mass. Thus in what follows we shall restrict to the study of the fundamental branch of solutions with 
a monotonic behaviour of the nA gauge function, i.e., w'(r) < everywhere. 

Some of the general features of the solutions are shown in Figure 4 where we plot the values Wh of the 
magnetic gauge potential at the horizon and Vq of the electric potential at infinity as a functions of rjj for 
a fixed k and several values of the charge parameter Q. For example, one can see that the solutions exist 



for w^ 1 < Wh < — 1, where the value w ™ increases with Q. For Q < one finds W ™ < 1, while 
w{ max) = 1 for Q > Q( c \ 



Q< c > = 16k. 
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Figure 3. The profiles of two EYMCS black hole solutions are presented as a function of the radial coordinate 
r. m(r) and cr(r) are metric functions, while V(r) and w(r) are electric and magnetic gauge potentials, respectively. 
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Figure 4. The value at the event horizon of the magnetic 
at infinity Vo (right) are shown as a functions of the event 
coupling constant n and several values of the electric charge 
of the electric charge parameter. 
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gauge potential w(rh) (left) and of the electric potential 
horizon radius rh for a fixed value of the Chern-Simons 
. Here and in Figures 5-8, = 16k is the critical value 



The behaviour of solutions as a function of the electric charge for a fixed event horizon radius (i.e., at 
fixed entropy) is displayed in Figure 5. The value of the CS coupling constant is also fixed there. One can 
notice there the existence of both a maximal and a minimal value for Q. The nA solution emerges as a 
perturbation of a RN black hole for a minimal value of Q given by riU(K). They exist up to a maximal 
value of Q = r\ + 16k where an extremal black hole is approached with w(r%) = 1. Also, above some value 
of Q close to Q^ c \ w(rh) becomes very close (but not equal) to 1. 

In Figure 6 we plot the order parameter J which enters the first relevant term in the large— r expansion 
of the magnetic gauge potential as a function of the Hawking temperature (i.e. a varying rh) and several 
values of the electric charge parameter. Again, the behaviour of J depends crucially on the value of Q. For 
Q < Q^ c \ J approaches a constant value as Th — > 0. The behaviour is different for Q > Q^ c \ J in that case 
increasing strongly with Th and taking always large values, which makes difficult its accurate computation. 
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Figure 5. Several relevant parameters are plotted as a function of the electric charge parameter Q for black hole so- 
lutions with a fixed value of the Chern-Simons coupling constant n and a given event horizon radius. The value of the 
mass parameter Mo and of the Hawking temperature Th are normalized with respect the critical Reissner-Nordstrom 
solution where a branch of non-Abelian solutions emerges as a perturbation. 
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Figure 6. The order parameter J which enters the asymptotics of the magnetic gauge potential at infinity is 
shown as a function of the Hawing temperature Th for a fixed value of the Chern-Simons coupling constant k and 
several values of the electric charge. Note that Th is normalized with respect the temperature of the critical Abelian 
solution where a branch of non-Abelian solutions emerges as a perturbation. 



4.3.2 The thermodynamics of solutions 

The nA solutions appear to exist for values of the CS coupling constant k > 1/8. Similar to the RN case, 
the nA black holes with given n, Q are found only for a finite interval of Th {i.e., of the entropy) only. The 
detailed picture depends however on the ratio Q/Q(°\ (We recall that for D — 5, = 16k.) 

For fixed Q ^ Q^ c \ the behaviour of the solutions is rather similar to the Abelian case and the temper- 
ature reaches its maximum at some intermediate value of the event horizon radius, an extremal black hole 
being approached for a minimal value of r/,. A plot of the horizon area as a function of the temperature 
reveals the existence of several branches of nA solutions. The typical picture can be summarised as follows. 
For a given k > 1/8 and any value of Q, a branch of non-Abelian solutions emerges as a perturbation of a 
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Figure 7. The scaled event horizon area an (left) and free energy / (right) are plotted vs. the scaled temper- 
ature tH for the non-Abelian solutions with several values of the electric charge and a given value of k. The branch 
of Reissner-Nordstrom (RN) solutions is also shown. 




critical RN configuration with = y Q/U(k). The entropy increases with temperature along this branch, 
which has, however, a small extension in both an and tu . This branch continues in a secondary one, where 
the temperature still increases, while the event horizon area decreses. Thus, these solutions have negative 
specific heat. For Q ^ Q^ c \ this branch ends for a maximal value of tu (whose value depends on Q) where 
a third branch of solutions emerge. This branch extends backward in (iff, Off) and has a positive specific 
heat. The Hawking temperature vanishes there for a minimal value r™ n ^ of the event horizon radius. 

As r —¥ , an extremal nA black hole solution with an AdS2 X S 3 near horizon geometry is approached. 

For Q < Q( c \ the parameters of this extremal black hole can be read from from (|3.15| . (|3.16|) . For example, 
the near horizon expansion of the solutions implies 

rt n) = | ~ wl) 2 , (4.6) 
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where Wh satisfies the cubic equation 

(64k 2 - w 2 h )Q 2 + 2/s(l + w h ) 2 (128K 2 (w h - 2) + w h (wl - 2w h + 3)) = 0. (4.7) 
For Q > Q( c \ the limiting extremal solution has 

r h = r ( h mm) = VQ-Q( C \ (4.8) 

and 

H(r) = -* (r - rh f + . . . , = ah (l + ^J^— X (r - r h r-^j +..., (4.9) 

2 

V(r) = (r — rji) + . . . , w(r) = 1 + wi(r - r h ) k + . . . , 

n L Oh 

where w\ and ah are parameters fixed by numerics and 

k = ~ (-1 + V5 + 32/t) > 1. (4.10) 

Although close to the horizon, due to the scaling relation (|3.6p . this extremal solution is essentially similar 
to the RN one, their bulk form is different. The nA solution presents a magnetic hair (outside the horizon) 
and a metric function cr(r) ^ 1. Another interesting feature of the nA configurations with Q > is 
that the magnetic flux lines are 'expelled' from the black holes as extremality is approached. That is, one 
finds that the nA magnetic field is vanishing on the horizon of the extremal black holes admitting the 
approximate expansion (I4.9[) . — 0. Thus these solutions seem to exhibit a sort of nA 'Meissner effect' 
which is characteristic of superconductiong media. This should be contrasted with the Q < case, which 
possesses a nontrivial nA magnetic field on the horizon. 

Some features of the nA black holes are shown in Figures 7 and 8 (left) where we plot the reduced area 
of the horizon an = 2TT 2 rf i /Q 3 / 2 as a function of the dimensionless temperature tu for a fixed CS coupling 
constant and several values of the charge parameter Q. The branch of RN solutions as given by (|3.38[) . f|3.41[) 
is also shown there. 

Furthermore, it turns out that the free energy F = A4 — ThS of a RN solution is larger than the free 
energy of a lower branch nA solution with the same temperature and electric charge, except for configurations 
with k close to 1/8 and small enough values of the charge, Q < /3. Therefore the nA black holes are 
generically preferred. These aspects are exhibited in Figures 7, 8 (right) where the dimensionless free energy 
/ is plotted as a function of the dimensionless temperature tu ■ Moreover, for the same values of the mass and 
electric charge, the RN solution typically has a smaller event horizon radius (and thus a smaller entropy), 
than the nA black hole [28]. Note, however, that most of the nA configurations have no RN counterparts, 
see the Figures 7, 8. 

Also, one can see that the interval of the scaled temperature tu covered by the set of nA solutions shrinks 
to zero as K approaches the minimal value 1/8. As k — > 1/8, all three branches described above collapse 
to a single point, which is the extremal black hole solution. This limiting solution admits a closed form 
expression and will be discussed separately. 

The overall picture is somehow different for Q = Q^ c \ in which case, despite the presence of an electric 
charge, the nA black holes behave in a similar way to the vacuum Schwarzschild-Tangherlini solution, with 
a single branch of thermally unstable configurations, see Figure 7. These solutions emerge again from a 
critical RN solution with = yj /U(n) and can be continued for an arbitrarily small value of the event 
horizon radius. As — > 0, the black holes with Q = approach a set of globally regular particle-like 
solutions, with tn diverging in that limit, as expected. 

4.4 On the existence of D = 5 stable black hole solutions 

Typically, the existence of an unstable mode of a nA configuration is associated with the zeros of the mag- 
netic gauge potential w(r). Thus the fact that we have found nodeless solutions suggests the existence 
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Figure 9. The potentials for the perturbation equation (|4.14|) are shown together with the corresponding mag- 
netic gauge potentials of the unperturbed solutions for several values of and fixed values of k, Q. For r% = 3.29, 
there is a negative region of the potential, V(r) < 0. 



of configurations which are stable against spherically symmetric perturbations. Moreover, the nA solutions 
are arbitrarily close to the Abelian RN configuration, which is known to be stable. 

Thus in this subsection we address the issue of linear stability of the black holes discussed above. For 
that purpose, we have to study the evolution of linear perturbations around the equilibrium configuration. 
For the pure EYM case and with a gauge group 577(2), this has been studied by Okuyama and Maeda 
in |21) . who derived the corresponding pulsation equation. Below we repeat their derivation with a slight 
modification due to the presence of a CS term and hence also an electric potential. Also, the purpose here 
is to show the existence of stable black hole solutions rather than study the unstable modes, which for our 
purposes here is of secondary importance. 

In examining such time-dependent fluctuations, we consider the following metric Ansatz generalizing 

ds 2 = — — - + r 2 dn\- N(r,i)a 2 (r,t)dt 2 , with N(r,t) = 1 ^-A (4.11) 

N(r,t) r l 

On the gauge field sector, we shall restrict our study to perturbations within the considered SO (4) x 50(2) 
model. For the U(l) part, one can take without any loss of generality an Ansatz with a single nonvanishing 
component V(r,t). 

The construction of a general time-dependent Ansatz for the 50(3) gauge group in D = 5 dimensions 
has been extensively discussed in |21) . Interestingly, it turns out that the corresponding YM Ansatz is much 
more restricted in this case than in D = 4 dimensions, containing only one potential w(r,t). This agrees 
with the result found by taking an 50(4) x 50(2) truncation of the general Ansatz (|2.14[) . (|2.15j) . 

Then the perturbed variables can be written as 

m{r,t) = m(r) + emi(r)e ,ffi + a(r,t) = <r(r)(l + eai{r)e tnt ) + (4.12) 
w(r, t) =w(r) + ewi {r)e mt + ..., V(r, t) = V(r) + eV x {r)e mt + ..., 

with m(r), <r(r), w(r) and V(r) a static solution and e a small parameter. 

After replacing in the general EYMCS equations (|2.7I) . one finds the following relations valid to first 
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Figure 10. A number of parameters of the D = 5 globally regular solutions are plotted as a function of the 
Chern-Simons coupling constant k. 
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(4.13) 



Thus the functions mi, Pi and o\ can be eliminated in favor of W\(r), leading to a single Schrodinger 
equation for w\: 



-y-o + UnX = O X 
art 



(4.14) 



where \ = w iVr and a new radial coordinate is introduced via 
potential in (|4.14p is 
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(4.15) 



One can verify that for n — V — 0, the above relations reduce to those found in [21] for D — 5 EYM system 
(note that, however, the unperturbed solutions there have infinite mass). 

The potential above is regular in the entire range — oo < < oo. Near the event horizon, one finds 
Un — > 0; for large values of the potential is positive and bounded. Standard results from quantum 
mechanics [44] further imply that there are no negative eigenvalues for O 2 (and then no unstable modes) if 
the potential Uq is everywhere positive. 

Indeed, our numerical results show the existence of black hole solutions with a positive potential t/(r*) 
for all value of k > 1/8 we have considered. As expected, all stable solutions we could find in this way have 
no nodes of the magnetic gauge potential w(r) (see the Figure 9 for such configurations). Therefore at least 
some of our solutions are linearly stable. The picture is however, quite complicated and depends on the 
values of k, Q and of the event horizon radius. For example, for solutions with k = 0.3 and Q = O.ISTQ^, 
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all solutions between the critical extremal black hole configuration (with ~ 0.702) and rt — 0.95, have 
U{r+) > 0. However, for the same value of the k, all configurations with Q = 1.05Q( C > have f/(r*) < 0. 

At the same time, we cannot predict anything if U(r+) is not positive definite. In this case we have 
to solve numerically the eq. (14. 14|) as an eigenvalue problem. Very likely, the full picture is complicated 
and a systematic study would represents a very complex task. We note only that, by using a trial-function 
approach (see e.g. |45j ) we could prove that a number of EYMCS solutions with one node in w(r) are indeed 
unstable. 



4.5 The globally regular solutions 

All D — 5 globally regular solutions emerge as zero event horizon radius limit of the black hole branch with 
Q = 0/ c ' . In principle, a disconnected branch of solitons with k < 1/8 may exist; however, we could not find 
such configurations. Heuristically, this can be understood as follows. Since, from (|3.11|) . the electric charge 
is proportional to the CS coupling constant, the existence of a minimal value of k means that below that 
value the electrostatic repulsion is too small compared with other interactions for a bound state to exist. 

The properties of these solutions are uniquely fixed by the CS parameter re and are somehow different 
from other nA solitons in the literature. For example, their mass is an almost linear function of re, while 
the value at the origin of the metric function g tt is very close to —1, see Figure 10. The shooting parameter 
b = —w"(0)/2 and the value at infinity of the electric potential take also small values. It would be interesting 
to find an analytical understanding of these numerical results. A typical D — 5 globally regular solution is 
shown in Figure 14 (left). 

Concerning the stability of solitons, the formalism proposed above for black holes applies also in the 
zero event horizon case. However, at this stage we cannot say something precise on their stability, since 
all globally regular solutions we have investigated have a negative potential. Also, the fact that there are 
no nodeless solitons (since to(0) = 1 and w(oo) = —1) strongly suggests that all particle-like solutions are 
expected to be unstable. 

4.6 The case k = a/8: An exact solution 

In what follows, we find it interesting to restore the (G, e)-factors in the general expressions. For re = a/8, 
one finds the following exact solution of the EYMCS equations within the metric Ansatz (|2.9|) employing 
isotropic coordinate^, and a gauge group SO(4) x SO(2): 

w(r) = 7rlJ' vir) = l\fS fi ^ ^) = ^)< /iM = ^ = 7(7)> ( 4 - 16 ) 



with f(r) 



Q^ 2 ( Q J 2 



2r 2 V<9 (c) (2r 2 + jy 



where Q > and J are arbitrary parameters, and — 16re/e is the critical charge parameter, — 
This describes an extremal black hole with nA hair, the regular event horizon being at r = (in 
these isotropic coordinates). The mass, electric charge, chemical potential and entropy of this solution are 



eVG ~ " e 2 ' e ' 4 V ttG e0>) ' 



15 In principle, this solution can also be written in the Schwarzschild-type coordinate system II2.10J I, For example, the relation 
between the Schwarzschild radial coordinate f and the radial coordinate r in l|4.16|l is f = r/ \J f(r). However, this results in a 
much more complicated expression of the solution. 
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such that Ai = $Q. This exact solution has a number of interesting properties. For example, one can show 
that the horizon is regular, with the following behaviour as r — > 

AM = % = iQ 'V Q( " + 1 + ^ + AM = (Q _ q 4 ,i; )2qM2 + °M), (4-is, 

for the metric functions, and 



for the curvature and Kretschmann scalars. Also, the magnetic field is vanishing on the horizon (this property 
is also shared by the extremal black holes with Q > Q^> and k > a/8), with w(r) = 1 — 4r 2 /J + 0(r 4 ) . 
The RN solution is found by taking J = in (j4.16[) . The mass, electric charge and chemical potential of 

this solution are still given by f|4. 1 T[) . while the entropy of the RN solution is Srn — Ja{ ^2 ' ) 3 ^ 2 ; which 
is higher than the entropy of the corresponding nA configuration. Thus the extremal Abelian solution is 
thcrmodynamically favoured, a feature which seems to be shared by other extremal solutions with k > a/8 
we have constructed numerically. 

In the limit Q = Q^ c \ the solution (|4.16l) describes a particle-like soliton, in which case 



f(r) = 



(J + 2r 2 ) 5 



(4.20) 



while, from KWi . M = 87r 2 /e 2 , Q = 32 v / G7r 5 / 2 /( V / 3e 2 ), and $ = ^/3/(ttG)/4. One can show that r = is 
a regular origin, with the following behaviour in that limit 

fo 20^ 2 J 2 

Mr) = § = 1 + -^y- + 0(r 2 ), / (r) = (J + 2Q(e)2)2 + 0(r 2 ), (4.21) 



and 



c)2 



5760gW 4 



One should note that the parameter J which appears via the magnetic gauge potential w(r) does not enter 
any global quantity (although these nA solutions are supported by a nonzero J). However, J has a physical 
meaning, since it enters also the metric components. Then, although the mass and charge are the same, a 
local observer could distinguish between the nA solution (|4.16p and the corresponding RN solution via the 
motion of a test particle. 

The reason for the existence of the exact solution (14.161) can be understood by noticing that the solutions 
reported in this work admit an interesting connection with a model considered in [29] . This connection 
follows from the observation that, for static configurations, the action (|2.ip with an SO(A) x SO(2) gauge 
group reduces essentially to a Einstein- Yang-Mills-Maxwell model, with a Chern-Simons-type coupling term 
between the U(l) and nA fields. Basically, with this group contraction down from 50(6), the CS density 
(12. 5|) reduces to the hybrid CS density in [25]. Moreover, as noticed in that paper, this model with k — a/8 
corresponds to the coupling of the super- YM theory to the D — 5 supergravity [46] . 

It is shown in )29j that for this value of the CS coupling constant, any flat space self-dual solution of 
the D — 4 YM equation with a gauge group SO(4) can be uplifted to D = 5 and promoted to soliton 
solutions of the full model. 



16 In |29_ in fact, the YM field is takes its values in one or other chiral represenation of S l O(4), namely for self- and antisclf-dual 
SU(2) fields. Here, by contrast, our "magnetic" YM connection is fully 50(4) valued, leading to the appearance of £5,6 = 75 
in ( t4^3l . 
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A slight generalisation of this construction can be summarised as follows. Let us consider a D = 4 
seed configuration consisting in a geometry described by a Ricci-flat line element da 2 (parametrized by the 
x-coordinates) with Euclidean signature, and a SO (4) nA field satisfying the self-duality eauationj^l 



F i:i = i\/deto : 7 5 ey fe; F w , 



(4.23) 



in a metric background given by da 2 . 

Then this configuration can be uplifted to solutions of the D = 5 EYMCS model in this work, for an 
SO(4) x 0(2) gauge group. The five-dimensional line element is 

ds 2 = -L-da 2 ~f(x)dt 2 , (4.24) 
/(») 

while the purely electric SO(2) potential reads 

y^ = \{^m. (4.25) 

In the above relations, f(x) is a solution of the Poisson equation 

V 2 (i) = -^iTr{i^}, (4.26) 

where the operator V 2 is taken with respect to the four dimensional metric da 2 . The D — 5 YM gauge field 
is F = F lJ dx l A dx 3 + (dV A dt) . 

Then one can easily prove that the full set of equations (|2.7p are satisfied, provided n takes the special 
value, 

e / 7rG a , , „„. 

^2V-=8- (4 - 2?) 

In principle, this approach can be used to uplift to D — 5 all four dimensional self-dual solutions of the 
YM equations, including configurations displaying no symmetries and multi-center solutions. Note that only 
soliton solutions of the equation (|4.26[) were considered in Ref. [33]. Here, this has been extended to the 
construction of black hole solutions by adding to 1// an extra part which is a solution of the homogeneous 
equation V 2 (y) = with suitable boundary conditions. 

The simplest case is found by taking da 2 to be the four dimensional Euclidean space and F^ the BPST 
instanton [47] . This leads to the spherically symmetric configuration (|4.16[) discussed above. However, a 
variety of other physically interesting solutions can be obtained in a similar way (this includes configurations 
with a D — 4 non-flat base space, e.g. the Euclideanised Schwarzschild metric and the D = 4 YM instantons 
in [48], [49] ) . Moreover, one can generalise this framework by including rotation in the D = 5 metric Ansatz, 
in which case it is possible e.g., to find a hairy generalisation of the BMPV black hole [50]. Such solutions 
are found beyond the simple Ansatz in this work and will be reported elsewhere. 



5 On D > 5 solutions 

Our results confirm the existence of both black holes and soliton solutions of eqs. (|3.1[) for D = 7, 9 as well. 
Then we expect that the EYMCS model considered possesses asymptotically flat, spherically symmetric 
solutions with finite mass for any D = 2n + 1 > 5. The numerical methods employed in this case are similar 
to those discussed for D = 5, though as expected the numerical difficulties increase with D. 

17 Thus all known D = 4, 5(7(2) YM instantons can provide a solution of ||4,23| I, 
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Figure 11. A number of parameters of D — 7 black hole solutions are plotted as a function of the electric charge 
parameter Q for fixed values of (k, r^). Here and in Figure 12 Q^ c ' = — 64k/5 is the critical value of the electric 
charge parameter. 
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Figure 12. A number of parameters of D = 7 black hole solutions are shown for fixed values of k, Q and a 
varying r h . 



Determining the pattern of the D > 5 solutions in the parameter space represents a very complex task which 
is outside the scope of this paper. Instead, we analysed in detail a few particular classes of D — 7, 9 solutions 
which, hopefully, reflect at least some of the relevant properties of the general pattern. 

Not entirely surprisingly, it turns out that the case D = 5 has some special properties. For example we 
have seen already that for D > 5, unlike in D — 5, the branches of nA solutions do not end in RN black 
holes. Also, according to our description of the heuristic argument cancerning Derrick scaling at the end 
of Section 2.2, finite mass/energy solutions can be constructed in spacetime dimensions D > 7 even in the 
absence of the gravitational term in the Lagrangian, i.e. for a fixed Minkowski background. 

The reason is that only in D = 5 is it necessary to have the Einstein-Hilbert term in the Lagrangian 
to satisfy the (heuristic) Derrick scaling requirement. In that case, the Yang-Mills term scales as L~ 4 so 
that in four spacelike dimensions the scaling L~ 5 of the CS term is not balanced in the absence of grav- 
ity, the latter scaling as L~ 2 . Generally, the CS terms scales as L~( 2n+1 \ in 2n spacelike dimensions. 
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Figure 13. A number of parameters of D = 9 black hole solutions are plotted as a function of the value w(rh) 
of the magnetic potential on the horizon for fixed values of ft, r^. 
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Figure 14. The profiles of typical globally regular, particle-like EYMCS solutions are presented as function of 
the radial coordinate r for D — 5, 7 spacetime dimensions. 



Then the Derrick scaling requirement can be satisfied if the Yang-Mills term scales as L 4p , provided that 

2n+l>2n>4p, (5.1) 

i.e., when D > 4p + 1. The Yang-Mills terms that scales as L~ Ap are the p— th members of the YM 
hierarchy [51] 

41-^Tr{F(2^}. (5.2) 
Then the Lagrangian that satisfies the requirement (|5.1|) is 

^matter = ^YM K ^CS ' n ^ ' (^-^) 

This is possible to satisfy for n > 3, i.e. D > 7. Thus for D = 7 and D = 9, the only choice of YM term is 
the usual p = 1 member of the YM hierarchy (i.e. the one in ([2TT|) ) , while for D > 9 either one or both of 
p = 1 and p = 2 YM terms (|5.2j) are possible choices, etc. All properties of these solitons are fixed by the 
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value of the CS coupling constant, the relation (|3.31[) being valid also in that case. A detailed discussion of 
the nongravitating YMCS solutions, including an existence proof, will be presented elsewhere. 

Returning to gravitating configurations, one notes that, however, to some extent, the families of D > 5 
black holes resemble their five-dimensional counterparts. In particular, the profiles of the solutions look 
similar to those exhibited in Section 4 and hence will not be plotted here again. As in D = 5, no multinodc 
solutions of the magnetic gauge potential w were observed for D > 5. There are, however, marked qualitative 
differences in some properties of the gravitation solutions in D > 5, vs. those in D — 5. These are discussed 
below. 

Starting with the dependence of black hole properties on the electric charge parameter Q for a fixed 
event horizon radius r^, we plot in Figure 11 the results of the numerical integration for agiven value of the 
CS coupling constant n. It seems that, for any r^, the solutions exist for a finite rangq^fl of Q. (Thus, as 
expected, one cannot find black holes with an arbitrarily large electric charge.) 

The limiting behaviour of these solutions at the limit of the Q-interval is qualitatively different from the 
case D = 5. In terms of the value Wh of the magnetic gauge potential on the horizon, the solutions exist 
for Wh € (u>h,min, 1) where the minimal value Wh t min depends on (although always with Wh <m i„ 7^ — 1, 
i.e. the branch of nA solutions does not join the RN configuration). The numerical results suggest that in 
this case the limiting configurations consist of extremal black holes with a regular horizon, see the Figure 
11 (left). Another feature of the solutions which is worth mentioning is that for intermediate values of the 
horizon radius there exists a large region of the parameter Wh for which both V'(rh) and w'(rh) are very 
close to zero. Moreover, for some intermediate values of Q, we notice the existence of two different solutions 
with the same charge parameter. 

The behaviour of solutions as a function of the event horizon radius is shown in Figure 12 for several 
values of the charge parameter Q and a given n. In terms of the scaled event horizon radius an and scaled 
temperature tu, the sets of nA black holes with Q ^ interpolate between two extremal configurations 
and one can notice the existence of three branches of solutions. The first branch of nA solutions starts from 
an extremal black hole with w(rh) taking a minimal value. This solution extends up to a maximal value 
of r^, where a second branch of nonextremal solutions emerges, extending backwards in r^. For the same 
r/j, the mass of one of these solutions is larger than the corresponding mass on the first branch (when they 
both exist). Also, this secondary branch has a negative specific heat, while the entropy increases with the 
temperature for the first branch. Finally, a third branch of solutions emerges for a critical r/j, which has 
again a positive specific heat. This branch ends in an extremal configuration with a minimal r^. 

The picture is different for Q = Q^°\ in which case the third branch is absent and the second branch of 
solutions extends to r/, — > 0, see Figure 12 (left). The limit of a vanishing event horizon radius corresponds 
to a globally regular, particle-like configuration. (This feature is similar to the D — 5 result.) 

The D = 9 solutions we have studied possess a similar pattern. A number of results in this case are 
shown in Figure 13 as a function of the parameter w{rh). One can see that the limits of the domain of 
existence corresponds to extremal configurations with Th = 0. 

We did not consider the issue of stability of D > 5 solutions. In principle, this is a straightforward 
extension of the work in D = 5, the problem reducing again to a single Schrodinger equation. The fact that 
for any D we have found nodeless solutions suggests the existence in all dimensions of configurations which 
are stable against spherically symmetric perturbations. 

As mentioned already, similar to D = 5, one finds also a different class of solutions describing globally 
regular solitons with Q = Q^ c \ In Figure 14 (right) we show the profile of such a configuration in D = 7 
spacetime dimensions. One can see some differences between this solution and the one in D = 5, the 
distorsion of the spacetime geometry being more pronounced in the higher dimensional case. 

6 Conclusions and further remarks 

The main purpose of this work was to present a general study of static spherically symmetric, asymptotically 
Minkowskian, solutions in a simple EYMCS model in D = In + 1 dimensions. Our choice of gauge group 

18 Although the solutions in Figure 11 have Q > Q( c > , for other values of we could find solutions with Q < Q( c > as well. 
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is SO(D — 1) x 50(2), which is just sufficient to support a nonvanishing Chern-Simons density. While 
the smallest simplest gauge group with this property is SO(D + 1), the asymptotic analysis suggested that 
there were no well behaved solutions in that case. This is in contrast with the asymptotically AdS solutions 
encountered in [36], which were found for the full SO(D + 1) gauge group (with D = 5). 

The CS term allows to avoid the usual Derrick-type scaling argument against the existence of D > 4 static 
nA configurations with finite mass in the usual EYM system. This provides an alternative to the higher 
order curvature terms of the YM hierarchy employed previously to reach the same result |25] , presenting at 
the same time a richer pattern. 

While we could display a number of analytic results valid for any (odd) values of I? > 3, our main 
numerical analysis was restricted to the case of D = 5 and and to a lesser extent to D = 7. In addition, 
we could confirm the existence of solutions for D = 9 as well. One should emphasise that the properties of 
the EYMCS solutions in this work are strikingly different from solutions to other nA models without a CS 
term, considered in the literature. 

The main interesting results are for D = 5, in which case they can be summarised as follows: 

• The black hole solutions emerge as perturbations of the RN solution, which becomes unstable when 
embedded in a larger gauge group. 

• The nA solutions are generically thermodynamically favoured over the Abelian configurations. 

• Some of the nA configurations were shown to be stable against small perturbations. 

• A solution in closed form is found for the minimal value of the CS coupling constant, describing an 
extremal black hole with nA hair. In the absence of a horizon, this becomes the soliton found in [2"9"] . 

Our results in D = 7 and 9 differ somehow from those in D = 5. Most importantly, the RN solution in 
these dimensions does not become unstable when embedded in the corresponding (larger) nA gauge group 
(in fact this property holds for any D > 5). As a result, the nA black holes do not emerge as perturbations 
of the the RN solutions. Instead, for a fixed value of the electric charge Q ^ Q^ c \ they appear to interpolate 
between two nA extremal black holes. These qualitative differences are a consequence of the choice of EYMCS 
model made here, which in any case does not seem to be a consistent truncation of a known supergravity 
theory. However, it is possible to search for different EYMCS models in higher (than five) dimensions, whose 
solutions are likely to fulfill the properties of the D = 5 solutions itemised above. For this we note that in 
D = 5 the CS term scales as L~ 5 , versus the YM term, which scales as L~ 4 . In this respect, it would be 
useful if in the higher 2n + 1 dimensions where the CS term scales as L~( 2n+1 \ the corresponding YM term 
would scale as L~ 2n . This can be achieved only in D = 4p + 1 dimensions, where the CS term scales as 
L~( 4p+1 \ by replacing the usual YM term F{2) 2 (the p = 1 member of the YM hierarchy) with the p— th 
member of the YM hierarchy which scales as L~ Ap [37] . It is obvious that this can only be done when n = 2p, 
thus the properties itemised above cannot be duplicated in D — 4p + 3 dimensions. 

This remark applies equally to the last of the properties itemised above, namely that closed form solutions 
can be constructed in all 4p + 1 dimensions, for a special value of the CS coupling constant. While in the 
p = 1 case discussed here the (usual) p = 1 BPST instanton is employed, in the case of a generic p the 
instanton [5Tj of the p— th member of the YM hierarchy is employed. This results in a general class of p > 1 
exact solutions, whose properties are similar to those of the D = 5 configurations discussed in Section 4.6. 

It is obvious that the black hole solutions in this work violate the no hair conjecture, that is, two distinct 
solutions can exist for a given set of global charges. Moreover, some of the nA solutions are really clasically 
stable, because they have maximum entropy among the black holes with the same mass and charge. This 
behaviour is somehow similar to that found in |52j for a family of monopole black holes in D = 4 Einstein- 
Yang-Mills-Higgs system. There too, a branch of monopole nA black holes merges with the magnetic RN 
solutions, which is unstable for some range of the parameters. Thus, similar to the case in [S2], the hairy 
black holes in this work may be relevant for the issue of the final stage of an evaporating D = 2n + 1 RN 
black hole. 

In principle, the study in this work can be extended in various directions. For example, it will be inter- 
esting to consider more general asymptotics and solutions describing black strings and p— branes. Another 
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possible direction would be to include rotation. However, even for the case of static solutions approaching 
the Minkowski background at infinity, with a gauge group SO(D — 1) x SO(2), we expect to find a variety 
of interesting solutions. For example, we expect both EYMCS solitons and black holes to exist, which are 
static but not spherically symmetric. Indeed, such solutions were found in D = 4 EYM system |14j . 

We close our discusion with some comments on another intriguing feature of the D = 5 solutions. Despite 
the different asymptotic structure of spacetime and the different horizon topology, these solutions have some 
similarities with the colorful black holes with charge in AdS space [TB], [S3], [53], which provide a model 
of holographic superconductors. In both cases, an Abelian gauge symmetry is spontaneously broken near a 
black hole horizon with the appearance of a condensate of nA gauge fields there, leading to a phase transition. 
Also, one can notice a striking similarity of the J(Th) curves shown in Figure 6 with some of those exhibited 
in the literature on AdS holographic superconductors. 

It remains an interesting open problem to clarify if the asymptotically flat EYMCS black holes may also 
provide useful analogies to phenomena observed in condensed matter physics. The first step in this direction 
would be to compute the conductivity as a function of frequency. This is obtained by perturbing the YM 
fields around the horizon. However, given the presence of several branches, the general picture is more 
complicated for asymptotically flat solutions, already for the fundamental RN set of solutions. Also, in the 
absence of a cosmological constant, the gauge/gravity duality (which seems to provide the deep reason behind 
the connection between general relativity solutions and condensed matter physics) is not yet understood. At 
the same time, some of the features of the AdS holographic duals of superconductors may occur for other 
asymptotics as well, being generic properties of certain classes of hairy black holes. We hope to return to a 
study of these aspects in a separate paper. 
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